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AN INFRARED SURVEY OF NEUTRON-CAPTURE ELEMENTS IN 

PLANETARY NEBULAE 

N. C. Sterling 1 , and H. L. Dinerstein 1 
RESUMEN 

El resumen sera traducido al espanol por los editores. We present results from an ongoing survey of the 
infrared emission lines [KrIII] 2.199 and [SelV] 2.287 /mi in Galactic planetary nebulae (PNe). Krypton and 
selenium may be produced by slow neutron-capture nucleosynthesis (the "s-process" ) during the asymptotic 
giant branch (AGB) phase of PN progenitor stars, and brought to the surface by convective dredge-up before 
PN ejection. We detect emission from Se and Kr in 65 of 114 Galactic PNe, and use the line fluxes to derive 
ionic abundances. We employ ionization correction factors based on coincidences of ionization potentials to 
calculate total elemental abundances, and discuss photoionization models designed to test the veracity of these 
corrections. The derived Se and Kr abundances of our targets range from approximately solar to enriched by 
a factor of 5, which indicates varying degrees of dredge-up and s-process efficiencies in the progenitor stars. 
In PNe exhibiting emission from both Se and Kr, we find that the relative abundances of these elements are 
generally in agreement with predictions from theoretical models of s-process nucleosynthesis. We examine our 
results for correlations between s-process enrichments and other nebular properties, such as CNO abundances, 
morphology, and characteristics of the central stars. PNe with Wolf-Rayet central stars tend to exhibit more 
elevated Se and Kr abundances than other nebulae. Bipolar nebulae, believed to arise from the most massive 
of PN progenitors, may have lower n-capture abundances than elliptical PNe. 

ABSTRACT 

We present results from an ongoing survey of the infrared emission lines [KrIII] 2.199 and [SelV] 2.287 /im in 
Galactic planetary nebulae (PNe). Krypton and selenium may be produced by slow neutron-capture nucleosyn- 
thesis (the "s-process" ) during the asymptotic giant branch (AGB) phase of PN progenitor stars, and brought 
to the surface by convective dredge-up before PN ejection. We detect emission from Se and Kr in 65 of 114 
Galactic PNe, and use the line fluxes to derive ionic abundances. We employ ionization correction factors based 
on coincidences of ionization potentials to calculate total elemental abundances, and discuss photoionization 
models designed to test the veracity of these corrections. The derived Se and Kr abundances of our targets 
range from approximately solar to enriched by a factor of 5, which indicates varying degrees of dredge-up and 
s-process efficiencies in the progenitor stars. In PNe exhibiting emission from both Se and Kr, we find that 
the relative abundances of these elements are generally in agreement with predictions from theoretical models 
of s-process nucleosynthesis. We examine our results for correlations between s-process enrichments and other 
nebular properties, such as CNO abundances, morphology, and characteristics of the central stars. PNe with 
Wolf-Rayet central stars tend to exhibit more elevated Se and Kr abundances than other nebulae. Bipolar 
nebulae, believed to arise from the most massive of PN progenitors, may have lower n-capture abundances 
than elliptical PNe. 
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1. INTRODUCTION 

The neutron(n)-capture elements Se (Z = 34) 
and Kr [Z = 36) may be produced by the s- 
process in low- and intermediate-mass stars during 
the thermally-pulsing phase of the asymptotic giant 
branch (AGB). During this stage of evolution, free 
neutrons are released in the intershell region primar- 
ily by the reaction 13 C(a,n) 16 0, and to a lesser ex- 
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tent 22 Ne(a,n) 25 Mg. These neutrons are captured 
by iron peak elements, which are transformed to 
heavier species by a series of n-captures and j3- 
decays (Busso et al. 1999). The newly-produced s- 
process isotopes may be conveyed to the stellar sur- 
face by third dredge-up (TDU), in which the con- 
vective envelope plunges beneath the (inactive) H- 
burning shell after a thermal pulse, thereby dredg- 
ing up He burning products (e.g. 12 C) as well as the 
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s-process-rich material. TDU is therefore an impor- 
tant contributor to the Galactic chemical evolution 
of C and n-capture elements, but it is not expected to 
occur in stars with masses lower than ~ 1.5 M©, and 
may vary in efficiency in stars with otherwise similar 
characteristics (Busso et al. 1999). This mass cutoff 
implies that the majority of AGB stars will not un- 
dergo efficient TDU, since low-mass stars are favored 
by the initial mass function. 

Carbon and n-capture elements are expected to 
be enriched in planetary nebulae (PNe) if the progen- 
itor star experienced TDU. However, the abundances 
of these species in PNe are poorly known, as the C 
abundance is inherently difficult to measure in ion- 
ized nebulae (Kaler 1983; Rola & Stasihska 1994), 
and n-capture element abundances had been mea- 
sured in only a few PNe before our survey. In fact, 
due to the low abundances of trans-iron elements 
(< 1CP 9 that of H) and the consequent weakness of 
their spectral features, it was not until 1994 that n- 
capture element lines were identified in the spectrum 
of a PN (Pequignot & Baluteau 1994). This led Din- 
erstein (2001) to realize that two previously anony- 
mous features at 2.199 and 2.287 /im in the spectra 
of several PNe are in fact fine structure transitions of 
[KrIIl] and [SelV]. Sterling et al. (2002) and Ster- 
ling & Dinerstein (2003) detected Ge in absorption 
in the UV spectra of six PNe, and derived the Ge 
abundance in five of these objects, but this only ex- 
panded the number of PNe with measured n-capture 
element abundances to seven. It is not possible to 
determine the history of convective dredge-up and 
AGB nucleosynthesis in PN progenitor stars, or how 
these evolutionary traits are related to other observ- 
able characteristics of PNe, from such a small sample 
of objects. 

We have conducted the first large-scale survey 
of n-capture element abundances in Galactic PNe, 
utilizing the two near-infrared emission lines iden- 
tified by Dinerstein (2001). Our results show that 
[KrIIl] 2.199 and [SelV] 2.287 fim are in fact de- 
tectable in a significant fraction (~60%) of Galactic 
PNe. Since Se and Kr are not expected to be de- 
pleted into dust (Cardelli et al. 1993; Cartledge et 
al. 2003), the gaseous abundances represent the true 
abundances of these elements. Therefore, [KrIIl] 
2.199 and [SelV] 2.287 [im are valuable tracers of 
heavy element abundances and AGB nucleosynthe- 
sis in PN progenitor stars. 

2. OBSERVATIONS 

We have observed approximately 100 PNe in the 
K band with the CoolSpec spectrometer (Lester et 



al. 2000) on the 2.7-m Harlan J. Smith Telescope at 
McDonald Observatory. Each PN is observed with 
a 2"7 slit at a resolution of ~ 500, which is sufficient 
to resolve [Kr III] and [Se IV] from most nearby fea- 
tures. We expand our sample by including K band 
observations of about 15 objects from the literature 
(Geballe et al. 1991; Hora et al. 1999; Lumsden et 
al. 2001). In our full sample of 114 PNe, we have 
detected Se and/or Kr in 65, for a detection rate of 
nearly 60%. 

The only important contaminants to the [Kr III] 
and [Se IV] lines are vibrationally-excited lines of H 2 , 
namely H 2 3-2 S(3) 2.201 and H 2 3-2 S(2) 2.287 /an. 
In the 15-20% of our targets which exhibit H 2 emis- 
sion, we employ a high-resolution (R ~ 4400) setting 
to separate [KrIIl] and H 2 , and use this result to 
estimate the relative contribution of H 2 to the unre- 
solvable blend at 2.287 /tan. In the rare case that the 
2.20 /im feature is too weak to be seen in the high- 
resolution setting, we assume the maximum contri- 
bution of H 2 to the blends, e.g. fluorescent excitation 
(Model 14 of Black & van Dishoeck 1987). 

3. ABUNDANCES AND CORRELATIONS 

We derive Kr ++ and Se 3+ ionic abundances (or 
upper limits) in all objects in our sample. A five- 
level atom is used to determine the level popula- 
tions of Kr ++ , while a two-level system is used for 
Sc 3+ . Transition probabilities and collision strengths 
are taken from Bicmont & Hansen (1986; 1987) and 
Schoning (1997) and Butler (2005, in preparation), 
respectively. We especially note the recently calcu- 
lated collision strength for [SelV] 2.287 /an (7.71 
at 10 4 K) by K. Butler, which previously was un- 
known. This allows us for the first time to derive 
absolute (rather than relative) Se abundances, and 
therefore we can compare the relative enrichments 
of Se and Kr in objects which exhibit emission from 
both species. 

In order to determine the elemental Se and Kr 
abundances, it is necessary to correct for the pres- 
ence of unseen stages of ionization. We account 
for this by assuming that Se 3+ /Se « Ar ++ /Ar and 
Kr ++ /Kr w S ++ /S, based on the similar ionization 
potential ranges of these species. It should be noted 
that this is an approximate method, and does not 
take into account the detailed atomic characteristics 
of each ion, such as the photoionization cross section 
and the rates of various recombination processes. A 
more robust treatment of the atomic physics is pos- 
sible with the aid of photoionization models, which 
can be used to derive more reliable ionization correc- 
tion factors. We have updated the atomic databases 
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of the photoionization codes XSTAR (Kallman & 
Bautista 2001) and CLOUDY (Fcrland ct al. 1998) 
to include Se and Kr. We will use these codes to 
compute detailed models of about a dozen PNe (with 
various degrees of excitation) , and derive ionization 
corrections from the results. In this paper, however, 
we cite results using only the approximate correc- 
tions based on ionization potential energies. 

The derived Se and Kr abundances span a wide 
range, from slightly subsolar (by about 0.1-0.3 dex) 
to enriched by a factor of five. We interpret this 
range of abundances to be due to varying efficiencies 
of the s-process and TDU in the progenitor stars, 
with the lower end of the range probably a metallic- 
ity effect. We find that about 40% of the 65 targets 
which display Se and/or Kr emission arc enriched 
in these elements. Converting this number to the 
fraction of all Galactic PN progenitors which experi- 
enced the s-process and TDU will require a detailed 
study of the sample bias of our survey. 

In those objects which display both Se and Kr 
emission, Kr tends to be more enriched than Se (see 
Table 1). Theoretical models of s-process nucleosyn- 
thesis (Busso et al. 2001; Busso 2003, private com- 
munication; Gallino 2005, private communication) 
predict that in the expected metallicity range of our 
sample (-0.3-1.0 Z Q ), [Kr/Sc] = 0.0-0.5. The ex- 
act value depends primarily on the size of the 13 C 
pocket (i.e. the mass of the layer in which free neu- 
trons arc produced), which is a free parameter in 
the models, and to a lesser extent on the progenitor 
mass (with more massive progenitors having larger 
[Kr/Se]). The observed [Kr/Se] ratios in our targets 
are generally in good agreement with the theoret- 
ical predictions, which indicates that the approxi- 
mate ionization corrections we use in this paper are 
not likely to be wildly incorrect. 

We have searched for correlations between the 
derived Se and Kr abundances and other nebular 
properties, such as CNO abundances, morphology, 
and central star type. In Table 1, we display av- 
erage Se and Kr abundances for a number of PN 
subclasses, including those with H-deficicnt, C-rich 
central stars (Wolf Rayet type, or [WC], and weak 
emission line stars, or WELS); PNe with large N and 
He enrichments (Peimbert Type I; Pcimbert 1978) 
which may indicate more massive progenitors (> 2.5- 
4.0 M Q ); and PNe of different morphological types. 
We also show the average abundances for the full 
sample, compared to the solar abundances compiled 
by Asplund ct al. (2005). 

We find that PNe with [WC] central stars tend to 
be more enriched in n-capture elements than other 



PNe. This is true for both Se and Kr, and is likely 
caused in part by the deep mixing and heavy mass- 
loss that these stars experience as they lose their 
H-rich envelope (e.g. Blocker 2001). The situation is 
less clear for PNe with WELS nuclei, as they seem 
to be more enriched in Se than PNe with H-rich cen- 
tral stars, but not Kr (however, the Kr results in 
this case are weakened by small- number statistics). 
Also, bipolar PNe, which may arise from more mas- 
sive progenitor stars, tend to be less enriched than 
ellipticals. This indicates that the s-process may be 
less efficient in more massive AGB stars. Such a 
trend is also seen in Kr for Type I vs. non-Type I 
PNe, but the difference in Se abundances is not sta- 
tistically significant. 

Interestingly, no relation is seen between the Se 
and Kr abundances and the C/O ratio. This comes 
as a surprise, since carbon is expected to be dredged 
up along with n-capture elements by TDU. Indeed, 
the abundances of n-capture elements have been 
found to scale with C/O in AGB (Smith & Lam- 
bert 1990; Abia et al. 2002) and post-AGB stars 
(Van Winckcl 2003). The lack of a detected cor- 
relation in PNe may be due to uncertainties in the 
derived Se and Kr abundances, as well as in the C 
abundances (which are usually not known to better 
than a factor of two in PNe). 

4. CONCLUSIONS 

We report preliminary results from the first large- 
scale survey of n-capture elements in Galactic PNe. 
We have detected emission lines of Se and/or Kr in 
65 of 114 PNe, from which we have derived the ele- 
mental Se and Kr abundances (or upper limits) for 
all targets in our sample. The recently calculated 
[SelV] 2.287 /an collision strength allows us to de- 
termine absolute Se abundances for the first time. 
We find a range of [Se/H] and [Kr/H], from slightly 
subsolar to enriched by a factor of five relative to so- 
lar. We interpret this as a variation in the efficiencies 
of the s-process and convective dredge- up in the pro- 
genitor stars. The abundances reported in this paper 
are rather uncertain, but will be made more robust 
as the ionization corrections are studied with the aid 
of the photoionization codes XSTAR and CLOUDY. 
The large number of Se and Kr detections will al- 
low us to examine AGB nucleosynthesis in PN pro- 
genitors, its efficiency in different subclasses of PNe, 
and consequently the Galactic chemical evolution of 
heavy elements. 

This work has been supported by NSF grants 
AST 97-31156 and AST 04-06809. 
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TABLE 1 



SE AND KR ABUNDANCES VS. NEBULAR PROPERTIES 



Property 


Derived 


Number of 


Derived 


Number of 




<Se/H> 


Se Detections 


<Kr/H> 


Kr Detections 


[WC] 


3.4±0.5 (-9) 


15 


6.4±2.9 (-9) 


3 


WELS 


2.8±0.7 (-9) 


9 


3.4±1.4 (-9) 


3 


Non-[WC]/WELS 


2.0±0.2 (-9) 


3G 


3.8±0.4 (-9) 


18 


Type I 


2.3±0.4 (-9) 


16 


4.1±1.2 (-9) 


5 


Non-Type I 


2.5±0.2 (-9) 


45 


6.8±1.3 (-9) 


19 



Bipolar 2.3±0.4 (-9) 

Elliptical 3.3±0.4 (-9) 

Round 2.1±0.6 (-9) 

Irregular a 2.5±0.6 (-9) 



11 2.6±0.4 (-9) 5 

21 4.2±0.7 (-9) 11 

3 ■•• 

6 1.7 (-9) 1 



Full Sample 2.4±0.2 (-9) 60 6.2±1.0 (-9) 24 

Solar b 2.14±0.20 (-9) ••• 1.91±0.36 (-9) 

a Includes Point- Symmetric nebulae 

b From Asplund et al. (2005) 
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